Super-exponential cutoff as a probe of annihilating dark matter 
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The spectrum of gamma rays from annihilating dark matter is expected to have distinguishable 
features if the strength of the background signal is not too large. Here we investigate the possibility of 
fitting the hypothetical spectra of gamma-rays from annihilating dark matter at the Galactic Center 
by a super-exponential template. By generating and analyzing mock data for several on-going and 
future experiments (HESS, CTA and DMA), we find that preference for a super-exponential template 
(over an exponential) might favour the presence of dark matter annihilating to t^t~ and jJ,'^ jjT 
annihilation channels. 



PACS numbers: 95.35.-|-d 95.85.Pw, 98.70.Rz 

The discovery of new physics often relies on our our 
ability to distinguish signal from background. Much of 
the time we have little prior knowledge about the 'more 
interesting' signal part of the data. Sometimes, however, 
the 'less interesting' background turns out to be no bet- 
ter known than the signal that we are seeking. When 
both signals have similar magnitudes, one could make 
use of the spectral differences between the signal and the 
background if they were sufficiently strong. 

It is necessary to invoke the concept of dark matter 
to explain many unrelated phenomena, such as the po- 
sitions of the peaks of the CMB power spectrum, abun- 
dances of light elements, cosmic structure formation and 
measurement of mass-to-light ratios at various scales. In 
a number of theories, dark matter would manifest itself 
by self-annihilating to Standard Model particles, with re- 
sulting galactic halo signatures that include high energy 
particles and photons [l|. 

The observed signal from any particular object or 
region of the Universe would inevitably contain back- 
grounds. The density of dark matter particles is deemed 
to be the highest at centres of dark matter halos. The 
strongest signal is therefore expected to be observed at 
the centre of our own Milky Way halo. 

The study of spectral features of astrophysical signals 
in the light of interpreting them as the spectrum of anni- 
hilating dark matter has recently attracted much atten- 
tion. The most distinct signature could be provided by a 
line feature [2tJ1i which may have been observed in the 
Galactic Center region at 130 GeV. However, whether 
the putative line feature is a real physical signal or has 
an instrumental origin is far from clear, cf. [5|. Less dra- 
matic, but possibly more realistic, spectral features that 
are hitherto unobserved are the so-called cut-off features 
Q , and box-shaped spectra features 0, Q . 



Here we study the spectral properties of a hypotheti- 
cal gamma-ray signal of dark matter annihilating in the 
Galactic Center region, and demonstrate that a simu- 
lated spectrum of dark matter annihilation to hard chan- 
nels (such as T'^T" and ^~^ii~) can be better fit by the 
so-called 'super-exponential' fitting function rather than 
by the commonly adopted exponential fit. 

As a template dark matter annihilation model where 
the signal is potentially strong, we adopt the dark matter 
density spike that forms around a central black hole in 
a dark matter halo or sub halo |9l-ll]|. This is visible as 
a point source, but similar arguments would apply to a 
diffuse dark matter concentration. 

The point source at the Galactic Center, associated 
with the central black hole Sgr A*, has been identified in 
radio, infrared and X-ray bands, as well as in gamma-rays 
by HESS [12[ . In the energy range of HESS it appears 
as a point source within the angular resolution of the 
instrument {Afl ~ 10~^). The HESS collaboration iden- 
tified a possible exponential cut-off in the spectrum of 
the Galactic Center source at about 15 TeV [1J|, which 
could also be explained by a broken power-law. The 
gamma-ray emission from source HESS J1745-290 and its 
high energy cutoff was explained by astrophysical mech- 
anisms in terms of acceleration and diffusion of cosmic 
rays [14l [l5| , and alternatively in terms of annihilating 
dark matter [16l - [l8l | . The next generation ground-based 
telescope CTA will provide superior spectral and angular 
resolution [l9[. For completeness in our analysis we also 
invoke the concept of the more futuristic experiment, the 
Dark Matter Array (DMA) ^. With such experiments 
in perspective, we argue below that improved spectral 
measurement of the putative TeV cut-off could discrimi- 
nate between the astrophysical and particle physics mod- 
els. 



In this paper first we describe the technique used to 
simulate the data and after that we present the resufts 
of fits of the simulated data by exponential and super- 
exponential functions. 

The mean number of events in the energy bin Ei , given 
the flux (/)(£') registered by an experiment with an effec- 
tive collecting area A{E) during time Ai, is estimated 
as 



N, = Ai- 



AE, 



(e' -E)^ 



dE / dE' A^Q[E')(l){E')e 



where 5E is the width of the energy bin, di is the vari- 
ance of energy resolution. For the sake of simplicity, we 
take the effective area Aqq{E) = Aqq{E = ITeV) to be 
constant. Once we have the mean number of events, we 
generate 100 realizations of the observations for a given 
bin and calculate the statistical error assuming a Poisson 
distribution. The values of energy resolution, typical ob- 
servation times, effective areas and systematic errors for 
HESS, CTA and DMA are summarized in Table HI 



We simulate the flux measured from the Galactic Cen- 
ter region as a sum of the dark matter annihilation spec- 
trum plus some unspecified power-law-like background 
(I){E) = 4>xiE) + 4'bg{E). While ground-based Cherenkov 
telescopes cannot distinguish gamma-rays from electrons 
and positrons, we assume that the major part of the elec- 
tron and positron background can be excluded from the 
Galactic Center region based on the spatial template. 

If dark matter is made of particles of mass AI^ that an- 
nihilate with a cross section (av) , the flux of gamma-rays 
from a dark matter halo with the density distribution de- 
scribed by p{l, ip, 0) is given by 



ME) 



(av) dN^ 
87rM2 llE 



p^{l,i;,(t))dldn, 



(1) 



dN^ 



where -r^ is the spectral density of gamma-rays per 
annihilation. It is standard to assume that the dark 
matter density in the vicinity of the Solar system is 
Pq{Rq = 8.5 kpc) = O.SGcVcm"'^. For demonstration 
purposes, we use a generalized NFW proflle with 7 = 1.42 
or equivalently take the flux to be boosted by a factor 
of 6 = 400 over the standard NFW density proflle with 
7 = 1.0 (consistent with [131) ■ 



(jj^iE) = 1.1 X 10"^ 



av 



3 X 10-26cni3s-i 



/100GeV\ dN 



M^ J dE 



cm ^s ^ 



(2) 
(3) 



We note that this proflle and/or boost is equivalent 
to adopting a density spike, for which the conventional 
assumption is 7' = 7/3 within ^0.1 pc. We assume 



that the background is given by a power law 4'bg{E) = 
Ai,gE~^^si with Fbg = 2.7 [21|. We assume smaller back- 
ground values for future experiments (see Table U for de- 
tails). 

An exponential cutoff is a template often used for de- 
scribing rapidly decreasing functions. The form 0(-E) ~ 
AE~^e~^^^^''"^'> has been used to describe the gamma- 
ray spectrum of the source at the Galactic Center [23 |. 
The continuum spectrum of annihilation of dark matter 
into photons for most annihilation channels has a pro- 
nounced cutoff at the mass of the dark matter particle. It 
is tempting to model such an abrupt cutoff with functions 
that decay more rapidly than an exponential function. A 
'super-exponential' function (j){E) = AE~^e~^^/^''^*^ is 
a simple generalization of an exponential function. 

We generated datasets that could be hypothetically 
observed by HESS, CTA and DMA for t+t- , p+fi', 
W^W" and 66 channels for dark matter with M^ = 4 
TeV using the PPPC4DMID code [H]. After a few tests 
we found that for our choices of experimental parameters, 
the T'^T~ and p^ pT annihilation channels benefit most 
from using the super-exponential spectral templates, and 
we decided to work specifically with these channels. Us- 
ing the technique of "scalable" windows, we binned the 
data in N\, ~ 20 logarithmic bins between Ea and i?b, for 
a range of e = Eb/ Ea- As found earlier, it is propitious 
to keep the geometric mean {E = \/EbEa) in a fraction 
of the dark matter mass E = e~°-2^M^ Q. 

For each e in the interval one to a hundred and each 
dataset, we found the best-fit parameters of the expo- 
nential and the super-exponential fits, and calculated 
the goodness-of-fit parameter x^- The average values of 
X^/dof and /3 are presented in the top and bottom rows 
of Fig. [1] respectively. 

Since the super-exponential function has one extra pa- 
rameter /? over the exponential function, we expect that 
X^ values for the former are at least as good as for the 
latter. Indeed, other than the region 1 < e < 2 for 
HESS, the x^/dof is signiflcantly better for the super- 
exponential fltting function than for the exponential one. 
For HESS, the flts with the super-exponential function 
prove to be good in the range of 1 < e < 20 for //+/Lt~ and 
T^T^ annihilation channel, while for CTA, the t~^t~ an- 
nihilation channel can be well fit with a super-exponential 
only in the range 1 < e < 5, and over a somewhat wider 



range 1 < e ^ 20 for the p n~ annihilation channel. 
However, the DMA flts with the super-exponential func- 
tion are not as good {x^/dof w 10). The reason for this is 
that with such overwhelming statistics and systematics, a 
simple super-exponential will not longer be a sufflciently 
good template to describe the desired signal. 

The bottom flgurcs show that the average exponent 13 
of the best flts is more or less a monotonically decreasing 
function of e. It is possible to claim that for e > 3 and all 
considered experiments, the /x"*"//" annihilation channel 
requires a greater /3 to flt the observed flux. 

After conflrming that for certain values of ApL the 
flux hypothetically observed by ground-based telescopes 



obs. time (hours) effective area (10 rn ) energy resolution syst. error background normalization Apl (GeV ^cm s ^) 



HESS 
CTA 
DMA 



50 
100 
5000 



0.18 
2.3 
23 



15% 
9% 

7% 



10% 

5% 
1% 



10" 

10" 



TABLE I: Benchmark parameters of current and future air Cherenkov telescopes observations of the Galactic Center. 
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FIG. 1: Top row: the averaged x^/dof of the best fits as a function oi E = \/EiE2. The values of x^/dof corresponding to 
the fits of the spectra by an exponential function are represented by dashed lines, and that of the super-exponential function 
by solid lines. The black lines correspond to the t^t~ annihilation channel, and the blue lines to the fj,^ fi~ annihilation 
channel. The normalization of the background varies as described in Table |ll Bottom row: The parameter /3 of best fits to 
the simulated datasets of the dark matter signal and background averaged over realizations as a function of E. The black 
solid lines correspond to t^t^ annihilation channel, the blue dashed lines - to fi^ fJ," annihilation channel. For all these plots 



E- 



'M^. From left to right : HESS, CTA, DMA. 



can be well fitted for e ~ 10, we studied the goodness of 
fitting as a function of Ap^ , the power of the hypothet- 
ical background (having fixed e = 10). In Fig. [2l one 
can see that as Apl grows, the exponential fits become 
better, while the super-exponential fits become worse. It 
is clear that if the background flux is comparable to the 
foreground flux around the cutoff energy (the mass of the 
dark matter), the relevance of sharp- featured templates is 
diminished. In Fig. [3J left inset, the average exponents (3 
are plotted as functions of Apl ■ We note that for better 
experiments (higher exposures, effective areas and lower 
systematic error) there will be correspondingly larger /3. 
The average /3 values of the best fits to the iX^ n~ anni- 
hilation channel arc systematically greater than those of 



the T^T annihilation channel. 

Examples of best fits to particular realization of fluxes 
observed from the Galactic Center by HESS and CTA are 
presented in Fig. [31 center and right plots respectively. 

We have investigated the use of super-exponential ver- 
sus exponential templates for fitting the signal composed 
of the flux of gamma-rays from dark matter annihilating 
at the Galactic Center and an unknown astrophysical 
background simulated by a power-law. We found that 
if dark matter annihilates to fJ.^ n" or t'^t^ channels 
than the super-exponential template provides good fits 
and in most cases much better fits than the exponential 
template. It allows us to assert that if, in general, the 
exponential fits to the gamma-ray flux from the Galactic 
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FIG. 2: The averaged x^/^of of the best fits as a function of Apl, the normalization of tlie background. The values of 
X /dof corresponding to the fits of the spectra by an exponential function are represented by dashed lines, and that of the 
super-exponential function - by solid lines. The black lines correspond to r'^r^ annihilation channel, the blue lines - to ^^ jjT 



annihilation channel, e = 10, .E - 



'M^. From left to right : HESS, CTA, DMA. 
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FIG. 3: The averaged /? of the best fits as a function of ApL, the normalization of the background (left). The values of 
P corresponding to HESS are presented by dotted lines, to CTA - by dashed lines, and to DMA - by solid lines. The blue 
lines (top three lines) correspond to fJ,'^ i-i~ annihilation channel and the black lines (bottom three lines) correspond to t^t~ 
annihilation channel. E — e~°'^^M-^, e — 10. The best spectral fits to t'^t~ and ^'^ fj,~ annihilation channel for HESS (middle) 
and CTA (right). The solid (super-exponential fit) and dashed (exponential fit) black lines (top) correspond to the fits to the 
t'^t~ annihilation channel and the solid (super-exponential fit) and dashed (exponential fit) blue lines (bottom) correspond to 
the fits to the n'^ n" annihilation channel. 



Center are significantly worse than the super-exponential 
ones, then the signal can be plausibly explained by dark 
matter annihilating to fJ-^ n^ or t'^t^ annihilation chan- 
nels. The inverse statement is even stronger - if there is 
no gain by using a super-exponential template over an 



exponential one, we can exclude n n 
lation channels. 
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